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Abstract 

We consider the duality between the four-dimensional S-matrix of planar maximally super- 
symmetric Yang-Mills theory and the expectation value of polygonal shaped Wilson loops in the 
same theory. We extend the duality to amplitudes with arbitrary hclicity states by introducing 
a suitable supersymmetric extension of the Wilson loop. We show that this object is determined 
by a host of recursion relations, which are valid at tree level and at loop level for a certain "loop 
integrand" defined within the Lagrangian insertion procedure. These recursion relations repro- 
duce the BCFW ones obeyed by tree-level scattering amplitudes, as well as their extension to 
loop integrands which appeared recently in the literature, establishing the duality to all orders 
in perturbation theory. Finally, we propose that a certain set of finite correlation functions can 
be used to compute all first derivatives of the logarithm of MHV amplitudes. 
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1 Introduction 



1.1 Dual conformal symmetry and Wilson loops 

Among all quantum field theories in four spacetime dimensions, maximally supersymmetric Yang- 
Mills (SYM) theory, in the 't Hooft limit of a large number of colors, is a most remarkable one. 
For one thing, it is the first and still only known example of an integrable theory in four spacetime 
dimensions. This integrability has been exploited successfully, so far, in the context of computing 
the spectrum of anomalous dimensions of the theory JlJ[2] . This theory is also singled out as being 
the prototypical example of the AdS/CFT duality |3|, which has produced many insights from 
strong coupling that have been instrumental in these developments. 

Integrability structures arise in the study of scattering amplitudes. By "integrability" we mean 
an infinite dimensional symmetry algebra acts on amplitudes. Most of these symmetries are hidden, 
non-local, symmetries that are not obvious from the usual presentation of the theory. The simplest 
example of such a symmetry can be seen from the expression of the one-loop correction to the 
four-point 1-loop amplitude. Ignoring regularization for a moment, this can be written, for a color- 
ordered partial amplitude, as [4] 

,1-loop _ .tree v - 2 *r f _A {Pi + P2? {P2 + PZ? 

4 4 9 C J i2^q 2 {<l + Pl?{q+Pl+P2?( q -p,?- [ ' 

Besides being very compact, this expression exhibits a very remarkable symmetry. This becomes 
manifest when one introduces the so-called dual coordinates Xi which are defined implicitly through 

Pi= x i~ x i-V ( 2 ) 
in terms of which the integral becomes simply 

d A x (x\ — x%) 2 (x2 — X4) 2 



(27r) 4 (x — xi) 2 (x — X2) 2 (x — xs) 2 (x — X4) 
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(3) 



The remarkable symmetry is the inversion Xi — > Xi/xf, which sends (xi — Xj) —> (xi — Xj) /x 2 xj: 
Equivalently, the integrand is conformally invariant. This "dual conformal symmetry" acts in 
momentum space and is therefore quite distinct from the familiar conformal symmetry of M = 4, 
whose arena is the original spacetime. 

Further evidence for this hidden symmetry accumulated in the course of higher-loop compu- 
tations, where it was repeatedly found that all integrals appearing in loop corrections shared this 
property unM- Of course, scattering amplitudes are infrared divergent, and regularization, for 
instance using dimensional regularization, will break the symmetry. So in the end most of the sym- 
metry is broken at loop level. Nevertheless, the breaking of the (bosonic) symmetry occurs in a very 
restricted and controlled fashion which leads to powerful constraints on the answer. Specifically, 
the breaking is entirely captured by the so-called Bern-Dixon-Smirnov (BDS) Ansatz Vn. Upon 
dividing by this Ansatz, any (finite) remainder is exactly conformal invariant [8j. In particular, the 
BDS Ansatz provides the full all- loop answer for 4- and 5-points amplitudes, because nontrivial 
conformal invariants do not exist in these cases. 
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In parallel development Alday and Maldacena [9~] studied scattering amplitudes at strong cou- 
pling using AdS/CFT. They confirmed the BDS Ansatz at 4- and 5-points but observed deviations 
for sufficiently large number of particles [To]; deviations at 6 points were also deduced from inves- 
tigation of Regge limits [Tl]. These violations are described by conformal invariant remainders, 
consistent with the symmetry analysis just mentioned. 

The strong coupling analysis led to a very geometrical understanding of the dual symmetry: 
a certain T-duality 12 takes planar Af = 4 to itself, loosely speaking interchanging the dual 
coordinates and space-time coordinates. It converts a scattering process to the expectation value 
of a polygonal Wilson loop with vertices at the points Xj, 

M n (pi,...,p n )<3>W n (xi,...,x n ). (4) 



Successive segments of the polygonal contour have length pf; we refer the reader to 13 for an 
excellent review. The "hidden" symmetries of the scattering amplitudes have now become ordinary 
conformal transformations of this Wilson loop. 

It makes sense to compute the same Wilson loop at weak coupling. The result was found, by 
explicit computations, to agree with scattering amplitudes, but with a certain very specific set of 
helicity configurations, the so-called MHV amplitudes [l4j[T5] . This Wilson loop/MHV amplitude 
duality has been confirmed up to two loops for the six-point function |16f[18| . 

In parallel development, it was understood how dual (super)conformal generators act on the 
polarization data that defines a scattering amplitude, allowing to define the symmetry away from 
the MHV limit. Tree amplitudes in that case have been shown to be invariant under a full dual 



superconformal 19 -23 , while loops carry only a (mildly broken) dual conformal invariance 19,24 

A certain small piece had thus been missing for some time in order to bring this story to its 
full completion: a geometrical understanding of the symmetry for non-MHV amplitudes, similar to 
that afforded by the Wilson loop in the MHV case. One of the aims of this paper is to provide a 
super- Wilson loop that fulfills just that purpose^ 

1.2 Loop integrands 

A significant extension of the scope of the duality was proposed recently in |26| . There it was pro- 
posed that the duality holds not only for final answers, but already for suitably defined integrands. 
The relevant "integrands" have been defined in two recent papers, in two very different contexts. 

In the context of Wilson loop computations, Eden, Korchemsky and Sokatchev [27| introduced 
a "loop integrand" by using the so-called method of Lagrangian insertion. In this method, loops 
are computed by taking repeated derivatives with respect to the coupling constant, which in SYM 
brings down powers of the Lagrangian density. At I loops these are to be integrated over I copies 
of space-time, the integrand being given as a tree-level correlation function of the Wilson loop with 
I Lagrangian insertions. This way the computation of a Wilson loop is turned into a set of d 4 x 
integrals. By adding suitable total derivatives to the action density, a "crural" Lagrangian insertion 
involving F af3 F af) instead of F^F^ can be used, which was found in |27| to be greatly convenient. 



1 While this manuscript was being prepared, another Wilson loop, very likely to be equivalent, has been derived 



from a very different perspective by Mason and Skinner 25 
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Independently, Arkani-Hamed and collaborators 128] including the present author introduced 
a loop integrand in the context of scattering amplitude computations. At the simplest level, this 
"integrand" is just the familiar momentum-space one, obtained after summing over all Feynman 
graphs and doing numerator algebra. Normally, such an integrand is not quite well-defined, because 
loop momenta are only dummy integration variables: it is unclear how to add different Feynman 
graphs under the same integration sign. One could always shift one diagram relative to the others. 
This ambiguity turns out to be absent in the planar limit, where region momenta Xi can be used 
to canonically add diagrams. This way a gauge-invariant integrand is obtained. The work |28j gave 
an efficient recursive procedure for computing it, starting from manifestly on-shell and physical tree 
amplitudes. 

The two integrands are known empirically to produce the same integrals, and it is therefore 
very natural to ask whether they are actually the same. Both integrals are chiral (non parity- 
invariant): in the first case due to a technical choice, in the second case due to one considering 
MHV amplitudes as opposed to its parity conjugate MHV. Thus they do not obviously disagree. 



Remarkably, explicit comparison in 26 has revealed that they precisely agree! In this paper we 
aim to provide an explanation for this fact. The explanation will, in fact, lead directly to a Wilson 
loop with polarizations, and to a general proof of the duality. 

1.3 Super- Wilson loop and recursion relations 

Our starting point is the singularities of the integrand defined by Lagrangian insertions in \27\ . 
This is a rational function of the insertion point y, with simple poles of the form l/(y— x«) 2 where y 
becomes null-separated from a corner of the Wilson loop. These singularities can be understood at 
the operator level: a singular propagator knocks out one factor of F a p from a Lagrangian insertion, 

F 2 ~ 1 \ a \ p F a P(y) + non-singular. (5) 

[y x i) 

Here A is some spinor associated with the separation (y — Xi) that is becoming null. The resulting 
correlation function can be made gauge-invariant in a natural way, by adding two segments to the 
Wilson loop joining Xi to y and back. Thus the residue at the pole is a Wilson loop with two more 
edges and F a p inserted at a corner: 



\ / 

In the scattering amplitude context, exactly the same residue is interpreted as the forward limit 




of a NMHV tree amplitude with two particles added 28,29]. Thus Wilson loops, with suitable 
insertions, know about the forward limit of NMHV tree amplitudes! 

As we will see, the duality is not limited to forward limits: with suitable insertions along the 
contour, general NMHV tree amplitudes are also reproduced. Furthermore, there is no reason to 
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stop at NMHV level. We will define a super- Wilson loop that knows about all N fe MHV amplitudes, 
and agrees with all MHV-stripped super-amplitude M n where 

An = <12>(23>...(nl) X Mn - (6) 

This new sweep of super- Wilson loops depends on the same momentum and polarization data 
as scattering amplitudes do. Importantly, within this enlarged sweep of objects, recursion relations 
become possible that relate objects with t Larangian insertions to other objects with fewer La- 
grangian insertions but more edges. These recursion relations will be seen to be identical to those 



found in 28 , establishing that the two sweeps of objects are the same, to all loops, for arbitrary 
polarizations, at the integrand level. 

This paper is organized as follows. In section [2] we introduce the coordinates that prove to 
be particularly well-suited for our description of scattering amplitudes, the so-called momentum 
twistor variables. In section [3] we begin the construction of the super- Wilson loop through a process 
of "reverse-engineering the answer." We complete the construction in section[4]by requiring that the 
object be supersymmetric, in a certain specific sense. Then in[5]we give a simple example confirming 
that the super- Wilson loop is indeed able to reproduce 1-loop NMHV amplitudes. The remainder of 
the paper is devoted to establishing the recursion relations, at tree-level[6]and for loop integrands [7} 
The paper ends with some conclusions, together with Appendices that give more details about the 
Wilson loop operator and various computations. In a final Appendix we outline, based on physical 
arguments, how divergences are expected to factorize within the Lagrangian insertion procedure. 



2 Lightning introduction to momentum twistors 

There are many useful ways to describe tree amplitudes in M = 4 SYM. Since we are interested 
in the dual Wilson loop picture, we should employ those variables that make the symmetries of 
the Wilson loop most manifest. Incidentally, these variables are also the ones which produce the 
most compact analytic expressions: These are the momentum (super-)twistor variables introduced 



by Hodges 23 



Momentum twistors can be used in any planar theory with massless particles. Their main 
achievement is to explicitly solve both the momentum conservation and mass-shell condition con- 
straints, allowing for a constrain-free description of the external data. The constrains can be solved 
in two successive steps. The first step was already implicit in the introduction: one can trivialize 
momentum conservation in a planar theory by introducing variables Xj such that 

Pi = x i~ x i-v ( 7 ) 

Planarity provides the required ordering of the particle labels. In the scattering amplitude picture 
the Xi are called region momenta, while in the dual theory they are the positions of the corners of 
the Wilson polygon. 

To solve the mass-shell constraints (xi — Xj_i) 2 = we proceed to assign a "momentum twistor" 
to each particle. This is a four-component object 

Z i = \ ..a ) = ( „ i\a J ■ ( 8 ) 



,,9 / \ t ■ Q A 

r 1 ! J \ 4Cf 
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where the spinor Xf is the one entering the standard decomposition pf a = XfXf, for pf = 0. A" is 
defined only projectively and so is Zf: Zf and UZf are to be identified, for any nonzero complex 
number 

The fact that this does solve the on-shell conditions will be explained shortly. As a first check 
one can count the degrees of freedom. We have 4n momentum twistor components, subject to n 
projective invariances and 4 translation invariances for the Xi, leaving (3n — 4) degrees of freedom. 
The original scattering data, on the other hand, has 4n variables subject to n + 4 constraints, giving 
also (3n — 4) degrees of freedom. Thus all constraints have become symmetries. 

The inverse change of variable can be constructed explicitly 23 1: 



The symbol (ab) is the standard antisymmetric 2-bracket, defined using the upper two components 
A" fe of Z a and Z\, (and more precisely defined in Appendix). Another useful equation is 

- . _ (i-i + (i i+i)^-! + (j+i i-vn4 , im 

while the Lorentz-invariant distance is simply 

(xi -x k ) 2 = (ii+l)(kk+l) x (ii+lkk+1) (11) 

where (abed) is the totally antisymmetric contraction ei 1 i 2 i 3 i 4 Z t a 1 Z l b 2 Z^ 3 Z^ . 

Geometrically, the point Xi in Minkowski space is identified with the line in momentum 

twistor space, that passes through the points Zi and Zi + \. The condition that Xi and Xj be null 
separated is that (ii+ljj+1) = 0. Geometrically, this happens if and only if the lines (ii+1) and 
(ii+1) intersect. 

We can now see why the mass-shell conditions are solved: the lines (ii—1) and (ii+1) automat- 
ically intersect, the intersection being just the point Zi. Thus pf = (xi — Xj_i) 2 = for any value 
of the momentum twistors. 

Momentum twistors make conformal properties manifest, even though conformal invariance is 
not a prerequisite for their use. The conformal group in four dimensions is SO (4,2) which is realized 
linearly on the momentum twistors. More precisely, momentum twistors transform in the spinor 
representation of (the double-cover) of SO (4,2). Recall that we are describing the dual conformal 
transformations, which are those acting simply in the X{ coordinates of the Wilson loop and not to 
be confused with the ordinary conformal transformations (acting on the original spacetime where 
the scattering process occurs). 

Because the antisymmetric tensor cabcd is invariant, any expression made out of 4-brackets is 
conformally invariant, while any 2-bracket breaks conformal invariance. 



To describe N = 4 polarization states one introduces the expansion for the on-shell states 30 



\g+) + fjf\^ iA ) + ^fjffjf^AB) + ^eABCDfjfvFvfli'F) + J$\9l)- (12) 
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(Notice our nonstandard use of tildes in this expansion. Untilded rj- variables will be denned shortly.) 
Capitalized indices range from 1 to 4 and carry the SU(4)^ R-symmetry charge, yielding two gluon 
helicities, four chiral fermions, four antichiral fermions and 6 scalars. 

Similarly to the Xi coordinates, one can introduce Grassman anticommuting coordinates Q% such 
that Xfrjf = 0f A — 6f_i- Momentum super-twistors are then defined as 

A? 

Z a =\ 1 = 1 I • (13) 




etxi 



The inverse map is 



k {i-li){ii+l) y ' 

Similar to the Xi, the 6i trivialize super-momentum conservation X^A^, which is part of the 
supersymmetry algebra. A consequence is that 9 and r\ only make sense on the support of suit- 
able super-momentum conserving delta functions. That is why, whenever one works with super- 
momentum twistors, one has to work with the MHV-stripped stripped amplitudes M n introduced 
above. Dual conformal invariance of the tree-level scattering amplitudes in N = 4 amounts to M n 
being expressible, using Zi and r\i variables, using only four-brackets. 

3 NMHV tree amplitudes from Wilson loops 

The (super) momentum twistors make tree amplitudes extremely well localized in the coordinates 
Xi\ this is the coordinate space of Wilson loop and so this is particularly relevant for us. For instance, 



consider the NMHV super-amplitude. It admits the very compact form 19,23 



M NMHV = £ [lii+ljj+l] (15) 

Ki<j<n-1 

where the square bracket is the ubiquitous supersymmetry invariant (Yangian invariant) 

\ ahcde \ = 8°^(r] a {bcde) + cyclic) 

(abcd)(bcde)(cdea)(deab)(eabc) 

With simple choices of fji the same expression gives NMHV tree amplitudes in pure Yang- 
Mills theory. The amplitudes with such values of rj are nontrivial linear combinations of those 
with prescribed rj values, and the conversion between r]- and fj- amplitudes requires full Af = 4 
amplitudes. In that sense the duality under consideration requires the full M = 4 supersymmetry: 
it is those amplitudes with prescribed rj that will appear on the Wilson loop side. 

To see how ^-components localize along the Wilson loop consider for instance the component 
proportional to f/iT^fc^fc+i in the scattering amplitude. It comes from only one ternj^J 

[1 2 3 A; k+1] ~ 77i?7 2 %% +1 x {12 ^ +l) ■ (17) 

2 To avoid cluttering the notation, here we omit the SU(4)p^ indices on the variables rjf. Whenever we write a 
product of four 77's, it is understood that all 77's have different i?-charge. Here, for instance, we could be extracting 
the component ^l^I^fc^fc+i- 
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Manifestly, this is completely localized to the corners x\ and x k of the Wilson loop: by inserting 
scalar field operators at corners of the Wilson loop 



4>AB{Xi) A B 

and connecting two corners with a scalar propagator this will be reproduced 



(12)(fcfc+l)(xi-x fc ) 2 



mmvkVk+i x 



(12kk+l)' 



(18) 



(19) 





This completely fixes the dependence of a putative super- Wilson loop on the product r/^j+i, if 
the tree NMHV amplitude is to be reproduced. Any other dependence on rjir]i + i could not reproduce 
this simple result. 

Other simple lessons gathered from the NMHV formula are that no dependence can occur on 
products such as r]ir)i+2 (because this would generate terms not in (Jl5|)), and that any operator 
inserted on the edge pi can depend only on the single momentum twistor r/j. 

To find out what these edge insertions are for scalar field components, consider the component 
^VkVk+i of the super-amplitude 



[12 3/c/c+l] ~ vimm+i x 



(13fcfc+l) 



(12kk+l)(2 3kk+l) 



(20) 



According to the above analysis we must reproduce this using a propagator linking the corner 
(kk+1) to the edge p2, or its endpoints. A natural Ansatz for the edge operator on segment 2, 
which will be motivated shortly, is 



1 \2vD4>AB(xi + tp 2 ) A B 
at r^—< % % • 



(21) 

io ( 2v ) 

The derivative is required to match with the dimensionality of the corner term and its form will 
be explained in a moment. For the moment this Ansatz contains an undetermined spinor v a . To 
determine it we compute the partial amplitude rj^rjkVk+i it produces: 



dt 



\2vd 



1 



\2v(x k - Xl) 



(2v)(kk+l) (xi + tp 2 - x k ) 2 {v2)(k k+l)(xi - x k ) 2 (x 2 - x k f 



(22) 



This conformal properties of this expression are not obvious. To make such properties manifest 
we rewrite the right-hand side using momentum twistor variables: 



(vkk+1 [1)(2 3]) 
(12fcfc+l)(2 3fcfc+l)(u2) 



(23) 
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where the square bracket means cyclic sum. It now becomes apparent that the choice v a = Xf 
produces almost the right answer: 

(lfcfc+1 [1)(23]) = (13fcfc+l) (31) 
(12fcjfe+l)(2 3ife*;+l)(12) (12fcfc+l)(2 3fcA:+l) (23 kk+l)(12) ' 1 ' 

where we have used the Schouten identity. The remainder is localized on the corner (23), and can 
be dealt with by a suitable corner term. 

Some comments are in order. First, the integral has produced a rational function of the external 
data (e.g., there is no logarithm). This is surprising because integration often produces logarithms. 
Here this can be traced to the chiral derivative acting on <fi. In fact, a rational answer will only 
be obtained for derivatives of the form \2aD aa or \2 a D aa . Investigation of the latter case reveals 
that it is incapable of producing the correct NMHV answer, leaving the above Ansatz as the only 
possibility. 

Therefore what we have found so far is that the coupling of the Wilson loop to scalar fields can 
be uniquely determined by "reverse-engineering" the NMHV formula: 

E i f « WlP t£ ( r + * ) *V + * Mxt) wk (*' " sw) ' (25) 

It is clear that by considering more and more polarizations we could continue with this "reverse- 
engineering" process and determine couplings to fermions and gauge boson. The tree-level NMHV 
analysis, however, will always be determined by free-theory two-point functions and is incapable 
of constraining non-linear couplings. To fully determine the super Wison-loop this way, we would 
have to consider N 2 MHV tree amplitudes, or NMHV loop amplitudes, and so on. Fortunately there 
is a better way to proceed. 



4 The full Wilson loop from supersymmetry constraints 

Tree-level scattering amplitudes are known to be invariant under dual super-conformal transforma- 
tions. They are in fact invariant under a whole Yangian algebra, but here we should focus only 
on those symmetries which the Wilson loop formulation is expected to make manifest. These con- 
tain the Poincare supersymmetries q, q of the Wilson loop which are expected to match certain 
elements of the SU(4|4) supersymmetries acting on momentum twistors. Specifically, the expected 
correspondence is obtained by matching the quantum numbers: 

92 - (26) 

(Notice that the R-charges carried by elementary fields are opposite in the T-dual theory.) The 
theory is also invariant under super-conformal generators s Aa and s% which will be discussed below. 
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It is useful to note that qA does not actually act on the original space-time scattering data: from 
that viewpoint qA is simply a redundancy in the 9 variables, much like translation invariance of 
the X{. Thus the super- Wilson loop has better be invariant under these transformations. However, 
contrary to translation symmetry, the bosonic part § dx^A^ is not trivially invariant under qA- this 
invariance is a very nontrivial constraint. 

To proceed we write the general couplings in the form 



p e Io dt£i{t) -pgjy 1 dt£ 2 (t) y 23 . . . p e Io dt£„(t)y^ 



nl 



(28) 



where the £i(t) = £j(xj_i + tp,j) are edge operators and the Vjj+i are vertex operators. In line with 
the above analysis, we will impose that £{ depends only on rji, and that Va+i depend only on r/j and 
T/i+i. Ultimately this is justified by the fact that a solution will be found with these constraints. 
The edge operator begins with the gauge connection £j = — Pi-A. The variation of this term 

is nonzero, as said above, q^Pi-A = \XfX i ^^. This forces us to add a fermion term such that 

d 



(<?2 + co£A Q A )^ = 0: 



1 

Ic 

This is unique if we assume that £i depends only on rj. 



£i = -Pi-A+^-X i ^ A r j f + ... (29) 



4.1 The Wilson loop in component form 

The procedure is systematic and by repeating it one obtains terms of higher and higher order in rj. 
With some algebra we find 



2£i = -XiXiA+ —X^AVi + 7T2 — 7 TT ^r\ — ^ ^ 

Cq ^Cq \l 1 1/ 



2 /„■„■ 1\ 'n '« 



1 XiXi-iXi-\D^ A bcd 1 AjAi_iAj_iA.;_iZ)F a b c d ron\ 

-&4 eABCD (ii-i? Vi * * ~ W£ ABCD <ii=i>s * * * * ' (30) 

together with the vertex terms 

V -II 4*AB ( a R A B \ , . 

Vh+i - 1 + ~ 2^ur Vi * ) + ■ ■ ■ (31) 

The somewhat lengthy expression for the vertex term is reported in Appendix. 

The overall factor of 2 in the edge term is just from dx^A^ — > dx ^ A A aa : for simplicity all terms 
are written as spinor products. In each term there is a unique way to contract the indices, taken 
conventionally to be with a lower spinor index on the left and upper index on the right. 

The constant cq has appeared because the two supersymmetries act on different superspaces: 
it is not a-priori obvious how to normalize one compared to the other. The correct value can be 
obtained by matching with the NMHV tree amplitude. With careful normalization (which was 
omitted in the previous section) we find 

C °=(l^r • (32) 
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Our construction and verification that the Wilson loop is invariant under the claimed q°i su- 
persymmetry has proceeded by brute force application of the supersymmetry generators given in 
Appendix. More precisely, the edges transform by total derivatives, which are cancelled by vertex 
variations. The precise form of the total derivative is given in Appendix. It is also important to use 
the equations of motion. 



4.2 Symmetries 

What about other symmetries? Consider the generator sa- It is easy to compute the lowest term 
in the edge variation: 

(s% + coM?;A)(-*iM + -X^Avf + •••) = -{\^)\^a - + • • • = + 0( V ). (33) 

dr)f c 

Thus the lowest term of the variation vanishes; here we have used the definition of Since the 
higher-order terms in rj are obtained by repeated action of qA, and since qA and sa commute, this 
can be used to show that the full operator is invariant. More precisely, saW is annihilated by qA, 
vanishes in lowest degree, and depends on edge i only through rji. But there is no such thing, so 
saW = to all orders in rj. 

The same argument shows that the operator is (formally) conformal invariant. In fact there is 
nothing to compute, since we already know that the Wilson loop is conformal invariant to order rp . 
So K aa W vanishes in lowest component, is annihilated by q (using that [q^if] oc sa and that sa 
is already known to be a symmetry), and depends on edge i only through rji'. this implies KW = 0. 
We say "formally" because, contrary to qA and sa, K aa contains bosonic derivative which make 
it more singular at short distances. Indeed it is well-known that K aa becomes anomalous at loop 
level |8j; qA and sa do not become anomalous, however. 

Finally, one might ask whether the operator is well defined, given that it contains null-separated 
fields. It appears that built-in cancellation mechanisms prevent ill-defined contractions of null- 
separated fields from occurring in it. For instance, the scalar field component (j){y) proportional to 



772% that is induced from edge 2 and vertex 2, as seen from (20), vanishes when <j> becomes null- 
separated from the endpoints of edge 2. Thus no singular contraction involving this component of 
the is possible. In other cases supersymmetric cancellations can be seen. Indeed a formal argument 
based on supersymmetry can be given to argue the absence of problematic contractions: by using 
12 of the supercharges Z a ^j to set three adjacent r/ A to zero, any potential power-divergent 
contraction between the two vertices will be prevented. Now of course one could do this anywhere 
around the loop, suggesting no such contraction occurs anywhere. Another line of argument is that 
the operator is essentially generated by repeated action of qA on the bosonic Wilson loop: since 
one never worries about the power-divergent contractions of (only about logarithimc divergent 
contributions at loop level), this leaves little reason to worry about such contractions for the super- 
extension of the operator. For these reasons we will discuss such issues any longer. 

In summary, we have defined a super- Wilson loop operator that is manifestly invariant under 
q\ and s\ and is also formally conformal invariant. Its scalar component agrees with that obtained 
in the preceding subsection. It is not invariant under the remaining supersymmetries, as will be 
discussed below. 
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4.3 Connection with the work by Mason and Skinner: superconnection 



In a remarkable recent paper, which appeared while this manuscript was being prepared, Mason 
and Skinner |25] introduced a Wilson loop in momentum twistor space, whose expectation value 
they argued (using MHV rules in momentum twistor space |31| ) reproduces all- loop scattering 
amplitudes in N = 4 SYM. The twistor transform of this operator, which they also discussed, gives 
a Wilson loop in the same four-dimensional space-time as the Wilson loop we are considering. This 
space-time Wilson loop they expressed as the integral over a super-connection, 



1 



-Tr 



Ve -IA 



(34) 



where 



A = 



-m AB e A e B 



l -e ABC DD^ A d B e c e D + ^e ABCD DF9 A 9 B 9 c 9 D + ... 



-d9 J 



1<Pab0 B + \e A BCD^ B C e D + l -e ABCD F9 B 9 c 9 D + 
2 6 b 



(35) 



Reference 25 explained the general method for constructing this operator and gave also explicit 



expressions up to scalar components. Here we wrote down a few extra higher-order terms by 
"guessing". (In each term shown there is a unique way to contract the indices. The notation is the 
same as above.) 

We believe that the shown terms are still not complete, though, since we also find that, if the 
object is to transform by super-gauge transformations under supersymmetry, it should contain other 
non-linear terms such as {4> A c, 4>bd}(9 2 ) AB {@ 2 ) CD ■ The full operator thus presumably contains 
term going all the way up to order 9 s . It would be useful in the future to understand how the 
twistor transform generates such terms. 

The two super- Wilson loops look at first sight very different. One contains corner terms, the 
other does not. But this is not an invariant statement: one can always re- write corner terms as 
suitable total derivates integrated along edges. It appears that the superconnection in |25| achieves 
this in a particularly clever way: for instance its rjiTji component gives 




XiXi-iDcj) ^ Xi\ i+ iD(j) (i-li+l)4> 
' 2<i 2—1) 2(ti+l) 



Xih-xDcj) + (i-li+l)<j> d 



2(*i-l) (i-l»)(tf+l) dt 



(36) 
(37) 



where we have used the Schouten identity. The first term is precisely our edge term while and the 
total derivative gives exactly our vertex term! We have also checked other scalar components. The 
agreement is even more striking for products like rjirji+i, which are pure vertex terms for us but 
which appear in Mason and Skinner's edges as pure total derivatives. Given the role played by 
supersymmetry in both constructions, we expect complete agreement for the other terms as well. 



A second minor difference is that the connection in 25 contains no factors of the coupling. 
Getting rid of factors of the coupling in (28) is actually easy, if one rescales ip — > c^tp, 4> — > Cq4> and 
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ip — > Cq'0- (This alters the reality condition (ip 01 )* = ip a , but within a Grassman path integral it 
is not clear whether this has any significance.) To get rid of the coupling that multiplies the field 



strength operator one can introduce an auxiliary field G a ^ , a la Chalmers-Siegel 32 and replace 
F/cq — > G inside the Wilson loop. After all this the Lagrangian density (76) becomes 



\l GapGaP ~ l^^tr B ] - ^ab,4>cd][4> AB A CD }) ■ (38) 



An important observation is that the first line is an action for self-dual N = 4 [32] . In particular, 
at order g° (tree- level) , the super- Wilson loop depends only on the self-dual sector. This is quite a 
surprise, and was also noted in [25] . 

4.4 The chiral better-half 

Tree amplitudes are invariant under both chiral and antichiral supersymmetries. The above makes 
manifest the q% and supersymmetries, but what about the other chiral halff]^] 
If one acts with the other half one finds for instance 



o rJXrt + coft 28m = \*{\iFe A {t) + -{4> BC , + c A^ A ) + . . . (39) 



In particular, this is nonzero. However, how big is this? 

If one uses the field rescaling just described, one trivially finds that this operator is explicitly 
of order A. That is, it contributes to 0(A) to any expectation value. This way the super- Wilson 
loop makes completely manifest the fact that tree amplitudes are invariant under a larger set of 
symmetries than loop amplitudes. Physically, this enhancement must clearly be related to the 
self-duality property which holds at zeroth order, although it is not yet completely clear how. 

The breaking of dual superconformal invariance in this formulation is thus explicit and unrelated 
to regularization. It exactly parallels what which was observed at the level of the answer, in the 



original paper on dual superconformal invariance 19 . As an operator, the super- Wilson loop is 



invariant under only a chiral half of the supersymmetries. 

5 A simple 1-loop NMHV example 

To illustrate the formalism, let us describe a particularly simple helicity component of the 1-loop 
NMHV amplitude proportional to ??i??2 r ?3 r /5- It is easy to very from the NMHV tree formula that 
this component vanishes at tree-level. This ensures that the one-loop correction is finite. In fact 
its value can be extracted from the general results for the NMHV ratio function that have been 
obtained in [l9j[33], in the scattering amplitude context, 

log u 2 log — (40) 



(1235) ° ux 



3 The author is indebted to E. Sokatchev for emphasizing this question. 
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where 



(1234)(4561) (2345)(5612) (3456)(6123) 

Ul ~ (1245) (3461)' U2 ~ (2356)(4512)' " 3 ~ (3461) (5623) ' 

In the super Wilson loop, two diagrams contribute to this helicity configuration: 



(41) 



6 



6 




3 



3 



Notice that the same topology contributes to amplitudes at various loop orders depending on which 
fields enter it: the shown diagrams contribute to 2-loop MHV amplitudes, 1-loop NMHV amplitudes, 
and tree N 2 MHV amplitudes, depending on how many powers of rj are extracted. 

We would like to check that these diagrams indeed reproduce the answer known from the scat- 
tering amplitude side. Notice we have a six- dimensional integral: 2 over the edge insertions and 4 
over the position of the vertex. We can perform these integrations in various orders. 

For instance we can do the edge integrations first. The fermion field sourced by edge 5 is 



and a similar formula gives the fermion fields sourced by edge 1. Inserting these result into the 
Yukawa vertex the first diagram gives for instance 



This has the form of a d A x loop integral as typically found in the scattering amplitude context. With 
the aim of making conformal invariance manifest we rewrite this as a momentum twistor integral. 
Using the notation of [28] the two diagrams just become^ 



Notice that the sum of the two diagrams is manifestly finite (all soft and hard collinear regions are 
killed by the numerator factor) and can be thus done directly in four dimensions. The integral can 
be done relatively easily using Feynman parameterization, but this will not discuss it further here, 
as we will now turn our focus to the more "Wilson-loop" form. Finite integrals of this type will be 
discussed at length in a coming publication |34| . 

4 The symbol (ABijkCMmn) = {Aijk){Blmn) - (A B). 




(42) 




(43) 




(44) 
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In the Wilson loop formulation it is possible to integrate over the vertex first. Using the three- 
point function 

(wmw * ( j*^-^ yr (45) 

the first diagram is thus written as 

dadb- vTttt " 7\ ~~ W7 ; — W/ — r (46) 

o (X4 + ap 5 - X2) Z {{1 - b)pi +P2) Z {X4 + ap 5 - x e - bp 1 ) z (23) 

Converted to momentum twistor variables the sum of the two diagrams becomes 

/245fi\ r dr A dr B (1235) [°° dr A dr c {l2M) 

{ 'Jo (5A1B) (5A23) (lB23) vcg { 'Jo (5Al2)(5A3C)(123C) ieg 1 ' 

where Zb = Z$ + tbZ§ and Zc = Z\ + tcZq. The r variables are related to the a, b variables 
through a = (^j+^g^) ■ The advantage of the r parametrization is that it removes all 2-brackets, 
thereby making conformal invariance manifest; the r-contours go from to infinity with constant 
argument such that e.g. Ta-t^ ^ s rea ^ anc ^ positive. 

Notice that divergence cancels between the two terms, but the individual (single-logarithmic) 
divergences as tb — > and tq — > have to be regulated using dimensional regularization. This 
amounts to (IB 23) — > (1B23)( ^^23) )~ £ ' ^ ne integrals are not very difficult and are carried out 
in Appendix; the result 

(1^5) l0gn2l0g S (48) 
agrees precisely with the one from the scattering amplitude side, confirming the duality at NMHV 
1-loop. 

This confirms the duality, but also gives a way to generate many very interesting relations 
between integrals of very different type. 



6 BCFW recursion relations: tree-level 

A powerful tool in the computation of scattering amplitudes is the so-called BCFW on-shell recursion 
relation. Similar recursion relations should thus apply on the Wilson loop side. We will now see 
that this is the case. In fact, the recursion relations are so powerful that this can be viewed as a 
proof of the duality. 

Actually, the first thing to check is that the Wilson loop at tree-level is a rational function. We 
have seen in our tree NMHV example above how the specific, chiral, form of derivatives on the 
operators ensures this. 

One general argument for rationality at order A goes as follows. The argument uses the antichi- 
ral supersymmetries q A and s A , which are valid at this order. They imply that the Wilson loop at 
this order is invariant under the full superconformal algebra (32 real supercharges); however "there 
is no way to supersymmetrize a logarithm": all such invariants are rational functions |35|. As will 
be noted below, it seems conceivable that a much weaker argument based on the fact that we are 
computing a tree-level correlation function in a conformal field theory would suffice, but this will 
not be pursued here. 
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6.1 Factorization limits 



Let us thus consider the singularities which arise when two corners of the Wilson loop, say x n and 
Xi, become null-separated from each other. Any singularity in this limit will be due to fields sourced 
at these corners, or on the edges very close to these corners. We can use supersymmetry to vastly 
reduce the number of fields which are sourced. Indeed, by using twelve of the Z a ^^ (which are 
true symmetries of the Wilson loop) we can put % = rji = rji = 0. We do not want to try to 
put r]i+i = as well, because the transformation which would achieve this becomes singular when 
(nlii+1) —7- 0, which is the limit we are interested in. However, there is no problem in putting 
three of the ?7's to zero. 

Let us be more explicit. Imagine we want to use supersymmetry to set r] n ,rji, rji and rjk to zero, 
where k is some other particle's label. Then we can solve explicitly for the shifted rj 

_ (n 1 i k)rg + cyclic 
T,i ~ (nlik) ' { 1 

These obey rj n = r)[ = = r]' k = 0, yet the Wilson loop is unchanged W(rf) = W{rf). The choice 
k = i + 1 would clearly be singular in the limit considered, but any other choice is perfectly fine. 

With the if variables the computation simplifies dramatically. Only is sourced near the 
corner x n of the Wilson loop, in Feynman gauge, 

A»{x) ~ £ log(x - x n f x ^ - P " ) . (50) 

87T 2 \p n -(x-X n ) py(x-X n )J 

We are ignoring terms that are nonsingular as (x — x n ) 2 — > 0. Unless this hits some derivative there 
will be no pole in l/(xj — x n ) 2 . We need either a field strength operator inserted at the corner Xi, 
or a derivative of the field strength integrated over edges close to x%. Since ^ = 0, using the form 



(28), the only term is from terms integrated along segment 



, 4 / dt\ i+1 \i\i\i g 2 (nl)(ii+l) 1 (i J) rf i+1 



Vi+1 J 2(i*+l)3 DF{Xt+tPt+l) ~ A^ct{nJ){Jl){iJ){Ji+l){x n - Xl ) 2 (ii+1) 4 (51) 

up to non-singular terms. The spinor Aj is the one associated with the null vector (xi — x n ). We'll 
simplify this expression in a second. 

In general there are interaction vertices in the bulk. We have to describe how they affect the 
single propagator contribution we have just computed. Scalar propagators near the light cone were 
analyzed in [36], where it was shown that the singular behavior simply gets dressed by a Wilson 
line. The same will be true here. We cannot directly use their result because we are considering 
a gauge field instead of a scalar field, however a simple extension will suffice. Actually we will not 
need to work very hard because at this stage we only need a tree-level result (in [36] the all-loop 
story was worked out). (Loops will be discussed below.) 

The argument goes as follows. First one notes that any interaction vertex with no derivative, 
inserted along the propagator joining x n to Xj, will remove the singularity. Thus one needs only 
consider vertices with a derivative acting on the singular propagators. These derivatives are essen- 



tially proportional to the (nearly) null vector [x\ — x n ). All this is as in 36 , here we only have to 
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make sure that the index structure of the gauge three-point vertex works out properly. Actually 
this is trivial, because (ccj— x n ) vanishes when dotted into the induced field (50), or into the field 
strength operator at the other end of the propagator. Thus (xi — x n ) must be dotted into the 
external gauge field, which reduces the vertex to a simple eikonal one. (The same argument rules 
out the gauge-scalar-scalar vertex.) We refer to | 36| for more details about the computation starting 
from this point. 



Thus the singular propagator (51) just gets dressed by an adjoint Wilson line. This means that 



the tree-level correlation function becomes the product of two Wilson loops, in the planar limit, 
since there are no interactions between the two sides. 




[n lii+1] x 




Taking residue 1 / (x r , 



Xi) 



i+1 ^ ^ i 

) 2 ) for prettier notation and expressing the prefactor 



(|51|) in terms of momentum twistors, which is done by simple algebraic manipulations, we thus 

(52) 



obtain the singular behavior: 

W n (l,...,n) -> [n lii+1] x W i+ i(l, . . . , i, J)W n _ m (J, i+1, . . . , n) 
where Zj is the momentum super-twistor 

Zj = (kii+1 [n)Ziy 



(53) 



(On the pole this is independent of Z&.) This is exactly how the tree amplitudes M n factorize! This 
has been derived when Vn = Vi = rfi = Vk = 0- However, being an equality between manifestly 
supersymmetric objects, it holds for arbitrary values of the r\{. 
In this expression 



[n lii+1] = <5((nlii+l}) 



<5°' 4 ((A; n 1 i)^,+i + (k 1 i i+l)r/ n + (ki i+1 n)rji + (k i+1 n l)r?i) 
(knli)(klii+l)(kii+ln)(ki+lnl) 



(54) 



is a supersymmetric version of a cut propagator. (This is also independent of Z/., as follows from 
the fact that, on the support of the bosonic delta function, all momentum twistors lie in the same 
3-dimensional subspace.) 

6.2 Tree-level BCFW recursion relations 

We have just seen that singularities of the super- Wilson loop factor the same way tree- level scattering 
amplitudes do. From this it is a short step to introduce the BCFW recursion relations 1 37, 38 1 (see 
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also (20) for the supersymmetric version) . Let us continue to discuss it here using momentum twistor 
variables. One considers the (super-)deformation of the external data 

Z n = Z n + zZ n — 1 (55) 

where z is a complex parameter. In the Wilson loop picture this corresponds to displacing the point 
x n but not the other points. Because the Wilson loop at leading order is a rational function of the 
external data, there is no problem in evaluating it at complex momenta. 
Then one writes 

W=l—W{z) (56) 



2iriz 

where the contour is a small circle around the origin. We want to deform the contour and express 
the integral in terms of the poles of W(z). The special values of z to consider are: 

• z = 00. In this limit x n —> x n -\ and the Wilson loop reduces to one with one less particle 
and Z n omitted. 

• z = — • These are the factorization channels quoted previously. 

• z = — (^-11) ■ At this point x n approaches spacetime infinity. The Wilson loop does not 
diverge in that limit, so there is no pole there. 

Summing up the residues one thus find 

W n (l,...,n) = W n ^(l,...,n-1) 

n—3 

+ ^[n-lnlii+l]W i+ i(l,...,i,Ji)W n - i+ i(Ji,i + l,...,n- I, hi) (57) 

i=2 

where the arguments are momentum super-twistors with 

hi = (lii+l[n)n-i\, Ji = (n-lnl [i) (58) 



This is the same recursion relation as obeyed by M n (see |39| for its manifestly dual conformal 
invariant formulation and [28] for its version in momentum twistor space similar to here). It is 
valid at zeroth order in the coupling, corresponding to tree-level in the scattering amplitude side. 
At this order, W4 = M4 = 1. Thus W n = M n for all n. The super- Wilson loop computes all tree 
amplitudes of the theory! 



7 Loops from chiral Lagrangian insertions 

We now discuss loop corrections to the Wilson loop. It will be very useful to use the Lagrangian 



insertion procedure used recently by Eden, Korchemsky and Sokatchev [26,27 in this context. As 
we will see, this allows the loop corrections to be matched on each side already at the level of the 
integrand. The expectation value of the Wilson loop can be written as the path integral 



/ 



[dAdipdfle'P 1 d4xC W. (59) 
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The leading contributions to the super- Wilson loop (dual to tree-level scattering amplitudes) are 
normalized to be of order g°. This is ensured by the proper choice for Co- Thus loop corrections can 
be obtained directly by taking derivatives with respect to the coupling constant, which will bring 
down powers of the Lagrangian. Actually, this requires some care, since the derivatives will also 
act on the explicit factors of the coupling in the expression in W. Also, one has to be careful with 
contact terms in the Lagrangian insertion. 

The Lagrangian insertion procedure, in the context of supersymmetric field theories, was first 



introduced in 40 . It was used to prove non-renormalization theorems [41] and also as an efficient 
tool for 2-loop computations |42 |. Contact terms play an important role in the consistency of the 
formalism [41 43 . 



7.1 Generalities 

Let us first consider a simple example in X(p 4 theory. The expectation value of some operator or 
product of operators is written as the path integral (overall normalization omitted) 

(0[<j)}) = j [DtyU-^+f-v^Oty]. (60) 

In this form one clearly has the expansion 



» Jv 



yi:—,yi 



£\" ^' 4! 4! 



(61) 



where the bracket is to be evaluated to zeroth order in the coupling. This is just the usual pertur- 
bation theory. Now let us rescale d> — > -A- to bring the Lagrangian into a form closer to what we 

v A 

have to deal with in a gauge theory. The same expectation is now written 



W]e^«VO[i], (62) 
A 



where 

V = ^V-^ 4 - (63) 

The reason for writing the kinetic term in this particular fashion will become clear shortly. 

The answer being the same as above we must still be able to compute loops by taking derivatives 
with respect to the coupling. But now one also gets terms from the explicit dependence on A of 
the operator, as well as from the kinetic term in the Lagrangian. The equivalence with ( |61| ) would 
appear to be lost. Consider for instance the first derivative of the two-point function (properly 
normalized so that the tree- level contribution is of order A ): 

d <j)(x\)4>{x2), 4>(X 1 )4>(X 2 ) , . f 4 (j)(xi)(j)(x2)C{z) 

'{ 7v 7v~) = -( ) -id y( ). (64) 



dX' VX V\ A 2 J A 3 



For a one-loop computation this is to be evaluated to order A . Now the first term would appear to 
be of order A -1 , however, the second term contains contact terms since ((D 2 (j)(y) — i0 3 (y))(/>(x)) = 
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i\5(x — y). These two sorts of terms are easily seen to exactly cancel each other. Furthermore, this 
use of the equations of motion flip the sign of the potential term, recovering precise agreement with 



the 1-loop term in (61) 



The same works at higher- loops. Derivatives of the explicit 1/A 2 factors from previous La- 
grangian insertions get cancelled by contact terms where multiple Lagrangian insertions go on top 
of each other. 

The moral is that, in this theory, by writing the kinetic term in the appropriate fashion it is 
possible to cancel contact terms against explicit coupling constant factors. Then it is correct to 
freely use the equations of motion inside the Lagrangian insertions. 

7.2 Chiral Lagrangian insertions 

All this carry over to any gauge theory, and so we return to M = 4 super- Yang-Mills for application. 



The Wilson loop ( 28 ) depends explicitly on the coupling constant through the factors c ~ A 1 / 4 . 
We would like to cancel this dependence against contact terms. This can be canonically achieved 
by writing 3/4 of the fermion kinetic term in the form (Dvp)^, and the remaining 1/4 in the form 



^(Dip). The precise formula is in Appendix (76). Similarly, using the self-dual form for the gauge 



kinetic term and the proper form for the scalar kinetic term given there, one finds that all contact 
terms in the Lagrangian cancel one to one against the explicit coupling dependence of W n . 
Thus the following expansion is valid 



(W n ) = "£(i\Y [ 



1 I Tit T ^OIl- shell {Vl) -Con-shell 
\ " n 



(.V n A 2 A 2 



(65) 



'yi,—,Vi 

where the bracket is to be evaluated to tree-level (A ). The on-shell Lagrangian is to be evaluated 



from (76) by freely using the equations of motion, dropping all contact terms. This gives 



£on-shell = \F aP F a ? - -^AB^tr 13 ] ~ ^ [<t>AB , C d] [<t> AB , <f°] ■ (66) 

Remarkably, this also happen to be a chiral operator: (j^vCon-sheii = 0, without any total derivative, 
as is easily verified]^] 

This is also obviously a conformal primary of dimension 4 (this especially obvious since, by virtue 
of the Lagrangian insertion procedure, it is only ever needed to tree- level order). Thus it is also 
invariant under the superconformal generators S°i. Under the remaining half of the supersymmetries 
q£ and it is not invariant, but transforms by total derivatives, for instance 

•^on-shell = -\D &a L$F<# + ^ Ba [4> AC , <f> BC ]\ • (67) 

We note that this is not quite the same as the Lagrangian insertion used in |27| or in |26|. It 



is however very simply related to, for instance, (3.16) in 36 . One only has to rescale the left 



handed and right-handed fermions separately in (3.12) and remove the anti self-dual part of the 



This chirality property also ensures that contact terms associated with coincident Lagrangian insertions produces 



terms proportional to £ on -sheii, guaranteeing the validity of the formula (651 to all loops 
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field strength by adding #-term. Complete equivalence of the results is therefore expected]^] We 
insist, however, that the form (65) is essentially forced upon us by the coupling constant dependence 
of the super Wilson loop. Any other writing of the kinetic terms would lead to uncanceled contact 
terms. 

These cancellations are suggestive that a field redefinition can remove the coupling dependence 
in W altogether, as in the above </> 4 example. In fact, this field redefinition is the Chalmers-Siegel 



one already mentioned around (38). We could have derived the above expansion quite directly 
starting from that form. 

7.3 BCFW for the loop integrand 

The first important remark is that the Lagrangian insertion procedure yields an integrand that 
is a rational function of its arguments yj. This was found in explicit computations [27]. In that 
context this was more or less guaranteed by the formulation as a tree-level correlation function in 
a conformal field theory. 

In our O(A ) analysis above, we used dual superconformal invariance to argue for rationality, 
but it appears a bit more difficult to make such an argument rigorous here since the Lagrangian 
insertion is invariant only up to a total derivative (although it is expected physically that this is 
always the total derivative of a rational function, as found in the scattering amplitude side (28). If 



this could be proved directly the argument would thus carry over.) In 27 , the main driving factor 
for rationality seemed to have been that the integrand was identified with a tree-level correlation 
function in a conformal field theory; thus it is also conceivable that a proof based on the weaker 
condition of conformal invariance is possible. 

This being said, it appears very likely that the integrand defined by Lagrangian insertion is 
a rational function in general. In what follows we will assume that this is the case, which will 
be justified a posteriori by the remarkable self- consistency of the emerging picture. So we let the 
integrand be a rational function of the insertion points and external data. Then we can continue 
using the BCFW deformation technology we used at tree-level to determine it. 

The contributions discussed in the previous section all go unscathed: they produce products 
of lower-point amplitudes, only, now one has to distribute the Lagrangian insertions in all possible 
ways between the two sides. However one expects new contributions too: terms where the deformed 
vertex x n becomes null-separated from a Lagrangian insertion point. In the scattering amplitude 
context [28] these terms are called single-cuts and they coincide with the forward limit of lower- loop, 



6 Our form seems also to have been used in 144], as the top component of the superfield used there. I thank 
E. Sokatchev for this observation. 
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higher-point scattering amplitudes. 




These "single-cut" terms were in fact our motivation for defining a super- Wilson loop, so that 
it could reproduce those terms. Let us thus conclude our story by verifying that the singularities 
of Lagrangian insertions indeed coincide with forward limits of W n +2 ■ 



7.3.1 Single-cuts 

This computation can again be greatly simplified if one judiciously uses supersymmetry to set 
rj n -i = 7] n = 7/1 = 0. In that case the fields which are singular at null separation from the corner 
(nl) are simply (50). Only the F 2 (y) term in a Lagrangian insertion will be singular, giving 

-F(y). (68) 



<? 2 1 (nl)X A X A 



16vr 2 x n -yf (nA)(Al) 

Interactions dress this propagator by an adjoint Wilson line, as verified above. Thus we have a 
fundamental Wilson loop with two more edges and a field strength inserted at the new corner, as 
claimed in the Introduction. 

In the BCFW formula the corresponding pole adds a new term which is a Wilson loop with a 
shifted point x n connected by an adjoint Wilson line to the insertion 

9 2 1 
167 2 " (ABn-ln)(ABnl) ABn—l l)(AB) 3 * X A X A F (y)- ( 69 ) 

Here the line (AB) corresponds to the point y and we have multiplied by a trivial factor 1/(AB) 4 
for future convenience (this factor enters the Jacobian J d 4 y — > cI^ZaZb/ (AB) 4 ) and 

n = (ABl[n-l)n), A = (n-ln 1 [A) B]. (70) 



This normalization for A is important here since (69) is not homogeneous in it. This computation 
is valid when rj n -i = r] n = r/i = 0. 

We have to compare this result with the forward limit of the super- Wilson loop ( 28 ) with n + 2 
particles. 

In the scattering amplitude context the forward limit is a somewhat singular limit but it is 
well-defined, for massless particles, at least in any theory with N = 1 supersymmetry 29 . A form 
particularly convenient for M = 4 was given recently in [281. This form uses a contour integral to 
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express the result in terms of strictly non-forward, hence manifestly well-defined, data. It looks 
even prettier when inserted into the BCFW formula, giving rise the BCFW term [28] 

/dVAs/ [A'B'n-lnl}M n+2 (l,...,h,A,B'). (71) 
J JGL(2) 

The contour integral is over GL(2) matrices G sending (A, B) to (A',B') = (A,B)G, and is to be 
done on residues around A' oc B' oc A. 



So one expects that the above expression ( 69 ) is equal to 



/AaAs/ [A'B'n-lnl}W n+ 2(l,...,h,A,B'). (72) 

J JGL(2) 

The verification is a straightforward computation using the explicit form of the edge and vertex 



terms. This computation is reported in Appendix with the result that expressions (69) and (72) 
agree precisely. 

7.3.2 Recursion relations 

This equality allows us to write the BCFW formula for the loop integrand defined by Lagrangian 
insertion in a form similar to that in |28| : 

W n (l,...,n;{AB} e ) = W n -x(l, n—1; {AB} t ) 

j n— 3 

+ J\ i+1 n ~ l n W i+^ •••>*>■£; {AB} L )W n - i+1 {Ji, i+1, . . . , n; {AB)} R ) 

ci i=2 



1 \^ 



/ [Aj Bjn+lnl] W n+2 (l,...,h AB] , Aj , Bj ; {AB}^ ) (73) 

JGL(2) 



/GL(2) 

where the shifted momentum (super- )twistors that enter are 

fij = {i i+1 1 [n— 1) n], Ji = (n— 1 n 1 [i) 
h AB] = {AjBjlln-ljn], Aj = (n-1 n 1 [Aj) Bj}. (74) 

In the second line the sum is over all ways of distributing the loop momenta {AB}i into two sets 
{AB} LjR . 



This is exactly recursion relation for the loop integrand found in 28 in the scattering ampli- 
tude context. It follows that, to all loop orders and for arbitrary external polarizations, the loop 
integrands for the super- Wilson loop and that for scattering amplitudes are identical. 

Both integrals are actually divergent. The Wilson loop is ultraviolet-divergent while scatter- 
ing amplitudes are infrared-divergent. The leading divergences on both side agree and are both 
controlled by the so-called "cusp anomalous dimension:" this quantity governs the singularities of 
null Wilson loops as well as the Sudakov factors in scattering amplitudes. On the other hand, 
subleading and finite parts need not agree, because natural regularization schemes on the two sides 
of the correspondence are different. To finish the proof of the correspondence, it is important to 
make sure that the logarithm of the amplitude (the relevant quantity) has a well-controlled scheme 
dependence. A general argument for this is presented in Appendix. 
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8 Conclusions 



We have constructed, directly from field theory, a family of supersymmetric Wilson loop operators 
in M = 4 Yang-Mills theory. This family depends on the same momentum and polarization data as 
scattering amplitudes do. By analyzing the factorization properties of this object in various channels 
at zeroth order in the coupling, we have derived recursion relations for it. These recursions are the 
same as those obeyed by tree-level scattering amplitudes, showing that the two objects are the 
same. Furthermore, similar recursion relations extend to appropriately defined "loop integrands" 
which can be defined on each side. The duality between the two objects thus extends to all loops, 
already at the level of integrands. 

The super- Wilson loop sheds some light on the nature of the duality. Understanding this at 
tree level is already nontrivial and interesting, since we have seen that in a sense loops are made 
by sewing trees. In the Wilson-loop side it turns out that only the self-dual sector contributes to 
tree amplitudes (this is especially clear in the Chalmers-Siegel form (38) and even more so in the 
work 1 25 1). This immediately makes a striking prediction on the scattering amplitude side: there 
must be a subsector of Af = 4 SYM that is sufficient to compute all tree amplitudes. This sector is 
certainly not directly the self-dual sector, because the asymptotic fields of negative helicity gluons 
are not self-dual by definition. However it should be somehow close to the self-dual sector. The 
amount of its non-self-duality is presumably related, through the T-duality, to the non-self-duality 
of the fields on the support of the Wilson loop. So at this stage it is not clear what the subsector 
is, and what its relation to self-dual J\f = 4 is, but it might be productive to understand this issue. 

It is important to note that we have not computed completely the form of the Wilson loop oper- 
ator: we have computed only those components that we needed, dropping certain nonlinear terms 
containing with 5 or more powers of 9 at the corners. Most of our results regarding factorization 
and recursion relations assumed little more than the existence of a supersymmetry-covariant Wilson 
loop, and supersymmetry always allowed us to make do without these terms. Nevertheless, since to 
our knowledge such a superconnection has never be written down in complete detail, we feel that 
it would be good to see it worked out fully somewhere. 

On a technical side, the examples that we have worked out have the unsatisfactory feature that 
conformal invariance is not manifest at all intermediate steps. The six-dimensional formalism for 
four- dimensional conformal field theories, presented recently by Weinberg |45| , might be very useful 
in that respect. Especially when combined with momentum twistors, which are in the spinor repre- 
sentation of SO(4,2) and allow for particularly compact expressions for chiral conformal invariants, 
this should provide a particularly efficient formalism. 

Another question on which the super- Wilson loop sheds some light is the longstanding obser- 
vation that tree amplitudes are invariant under a full Yangian symmetry, while loop amplitudes 
(even those which are finite) do not exhibit it [19] . In particular, half the dual supersymmetries 
d/dO are preserved to all loops, while the other half 9d/dx is broken starting from 1-loop. This 
has been understood as a holomorphic anomaly in [24]. The Wilson loop defined here and in [25] 
makes the unbroken half of the supersymmetry fully manifest but explicitly breaks the other half, 
As discussed in section (4.4), the Wilson loop, constructed from various other requirements, turns 
out to violate explicitly these symmetries but by an amount explicitly proportional to g 2 . Thus the 
empirical observation of [19] are built-in into the operator. The symmetry enhancement that occurs 



24 



at tree-level must clearly be related to the self-dual sector, although precisely how is not completely 
clear yet. 

The explicit breaking suggests a novel way of computing what has appeared as an "anomaly" 



before. Indeed, we have found that the variation is given by (67) 



iA* d ,__,„, , £ ^ ((c iP A + F9 A + . . .) a W n 



- a^-MWn) = f dx &a j^-y^ , (75) 

where the (color-adjoint) parenthesis is to be inserted along the Wilson loop; a similar expression 
holds for the superconformal generators s A . As just remarked, the right-hand side is manifestly of 
order g 2 . Furthermore, at any loop order, the right-hand side can have at most single-logarithmic 
divergences. In fact, it appears very likely that these divergences are concentrated to the corners 
of the dx integral, so that the correlation function itself is finite to any loop order. If this formula 
is correct, it might prove a particularly powerful way of computing derivatives of the logarithm of 
MHV amplitudes at higher-loop levels; these derivatives, and the full answer, are known to exhibit 
remarkable simplicity and structure [46] . 

The author would like to thank Nima Arkani-Hamed, Lionel Mason and David Skinner for useful 
discussions and Emery Sokatchev and Gregory Korchemsky for very useful correspondence and for 
a careful reading of this draft. Hospitality is gratefully acknowledged from the Perimeter Institute 
in Waterloo, Canada, for a wonderful stay during which part of this work was conducted. The 
author gratefully acknowledges support from NSF grant PHY-0969448. 
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A Lagrangian and supersymmetry generators 

We now describe, mostly for the sake of fixing our conventions, the well-known J\f = 4 Lagrangian 
in a form which will be convenient for us: 



- F aF a/3 + 



+ \<t>ABD 
o 



2 iAB 



+^ AB [^r B ] + -^ AB [^aM\ + 1~^abAcd][<P AB A CD ]- (76) 



—o 1 f d 4 xC 

This enters the path integral through the factor e ^ . The self-dual field strength tensor is 
defined through 

F aa = e ap F <*P + e «/9 F a/3 ( 77 ) 

where F^ v = d^A u — d u A^ + [A^,A U ]. (We suppress all color indices; LA^,^]" = f^A^A^ and 
The supersymmetry transformations, under which C transforms by a total derivative, are 



ffXAM = e a ^ A , q%^ B = D a ^ AB , (78) 

q * A <p BC = r [B 6° A \ = 5 B A F^ + \e^ AC , <t> BC \. (79) 

together with their parity conjugates, which are simply 

q A *4>BC = ^ q[ b^ v q A Hi = 6%-FW + \e^[<p AC , </> BC ]. (81) 

Our writing of the kinetic terms in the Lagrangian density may appear somewhat unusual; the 
motivation for it is explained in the main text. 

Let us finally summarize the rest of our notations. We have p a/3 = p'VjJ where a^ a = (i, —ia) = 
~ a ita- ^he s P m or products are (vw) = v a v a and [vw] = VaV a where indices are lowered using 
v a = v^epa and Va = v^e^ with e the antisymmetric symbol with e\2 = = 1- Indices are raised 

using the inverse symbol with e 12 = e 12 = —1. The fields obey the reality conditions (f)* AB = <p BA 
and (^ Aa ) f = 

B Vertex terms 



In this section we report our computation of the vertex term in the super- Wilson loop (28); for 
completeness we reproduce the edge term also. 

c 2cq 

1 AjAj_iAj_iDV A bcd I AjAj_iAj_iAj_iL>F A B c d / r,\ 
>e A BCD ... tto Vi Vi Vi ~ -rrn e ABCD ... lX o Vi Vi Vi Vi , (82) 



3!cf °" (ii-l) 2 h 11 h 4!cf °" <it-l)3 
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Va+i - 1 + 2 Vi+iVi - » / ■ , „ A Vi Vi 



c§(i+l*> V + 2 (i-lt) 

1 i-li+l)((^-l)A m + (ii+l)Ai_i)^ 1 A i+ i^ 



' 3!c§ (i-li) 2 (ii+l) 2 111111 ' 2c 3 Q (ii+1) 2 

1 Xi4> , 1 (i-H+1) 2 AB 4 

- 2c3 (ii+ip^ 1 ^ 1 + 44! (z i+ i)2 ^ * 

1 (i-H+l)((ii-l) 2 A m A i+ i + (ii-l)(i7+l)A i+ iAi_i + {ii+l^Xj-iXj-^F 4 
4!c^ (ii-l) 3 (ii+l) 3 ^ 

+0(r H r H+1 r } 2 ) + 0{rf). (83) 

We found it difficult to find a simple organizing principle for the vertex terms. Fortunately, the 
omitted terms in the vertices will not be need for the analyses in this paper. These omitted terms 
are all proportional to at least one power of rji and one power of r/i+i, and we expect them to form 
a nontrivial series ranging all the way up to r]frjf +1 . 

This is "supersymmetric" (under half of the supersymmetries) in the sense that edges transform 
by a gauge transformation (q^ + CoX a -^^)£i = (dt — [£%, )Xf A , where 

V a — (± B 1 Xi-\1p B c D 1 Aj-lAj-lF B c d\ 

XiA = co^T) - 2^ eABCD lJ^l)^ * ~ 3^ 6ABCD (ii-l)* * * * ) ■ 

The vertices transform by the same "gauge transformation," + cqX 01 ^g)Va+i + X^Va+i — 
Va+iX^A = so that the Wilson loop (28) is invariant. 



C Integral in the 1-loop example 

In this Appendix we consider the integral 



OO 



,>rf , ^^(1235) r dr A dr c {l2M) 

{ 'Jo (5A1B) (5A23) (lB23) reg { ' J (5Al2)(5A3C)(123C) reg 1 ' 

which arose in our NMHV 1-loop example, with Za = Z4 + taZq, Zb = Z2 + tbZ\ and Zq = 
Z2 + tqZ4 and the contours as described in the main text. The (single- logarithmic) divergences as 
tb — > and tq — > are regulated using dimensional regularization, (1B23) —> (1B23)( ^^^23) 
The divergences cancel between the two terms. After doing the tb and re integrals one is thus left 
with a finite integral 



amsrN r dTA 1 r ^ (5^12)(5A34)(6123) \ 
{ 'Jo (5A12)(5A23> g V(5A61) (5423) (1234) J 



(85) 



To proceed further we rewrite this as a function of cross-ratios by rescaling ta — > a (5612) ' sucn that 
the integral becomes 

1 ~ l/u 2 r da U{l + a) (4512) 2 (3456)(6123) \ 

(1235) J (1 + o)(l + a/u 2 ) ° g \l + a/u 2 (1234)(2345)(4561)(5612) y ' [ 1 
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We now observe that the rational factor has a symmetry under a —> u 2 /a. Upon symmetrizing the 
integrand, the a-dependence inside the logarithm disappears. Thus the integrand does not produce 
dilogarithms but only products of logarithms: 

l-l/«2 [°° da / (4512) (3456) (6123) \ 1 u 3 

(1235) J (1 + a){l + a/u 2 ) ° g V (1234) (4561) (5623) ) ~ (1235) ° gU2 ° g ^" ( ' 

D Forward limits 

In this section we check that the contour integral version of the forward limit 

I d% A d 4 r] B I [A'B' n-1 n 1]W„ +2 (1, . . . , h,A,B'). 

J JGL(2) 



reproduces just (69) 



A simple and correct way to do the integral over the GL(2) is as follows. Since invariance under 
the diagonal component of the GL(2) is already manifest, the integral to be done is really over 
only the off-diagonal components. This can be done in two steps. First one integrates over the 
off-diagonal component A' = A + t\B. The only dependence on t\ is through the R- invariant, 
which readily gives 



dTl \A'Bn In 11 = &Hva{B n n-1 1) + cyclic) 

2ni J {ABn-ln){ABnl){ABn-ll){Bn-lnl 



To proceed further we restrict to r? n -i = rj n = n\ = as in section (7.3.1). In this case n % = due 
to the fermion <5-function: all the r\A-,r\B dependence of W is through Zb> = Zb + tau^Z^- Thus 
we are extracting the component of order rj B ,. We can use this to do perform the fermion integrals, 



'4 a 5° l4 (VA(Bnn-ll) + cyclic) , 4 , 

d r] A d 7] B — — -Tj U| {r)B + T 2 r) A ) = (B n-lnl 

\B' n—1 n l) z 



2 



(90) 



The T2 integral is to be done on a small circle around the special value where (B'n—lnl) = 0. 
Now all we got so far near this point is a double zero. However some terms in the Wilson loop 



(82) contain denominators l/(ii— 1) —> (B 1 A), which produce poles. The terms in the Wilson loop 
contain at most three such powers, which is exactly what we need to obtain a finite residue. Thus 
we need only keep those terms with three powers of l/(i—li). These are 



/ dr 2 (B 1 n-lnl) 2 (Al) f 1 f 1 ~ , , . . , , 

i — ^ ^t^t4 x — / dt\t\i\i\iDF(y + tp) + \i\iF(y+p) 

J 2m (AB n-ln) (AB n 1)(AB n-1 1)(B> A)* (B'l) \ 2 J aaaa ^ v> a a W p) 

(91) 



The edge integral in the parenthesis is a total derivative, which cancels precisely against the vertex 
term at one endpoint. There remains only the other endpoint, located at y. Finally the contour 
integral can be expressed as 

TOn 2 (AB n-1 n) (AB n 1) ABn-ll){AB) 3 X X a X a F ^)- ( 92 ) 
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This agrees precisely with (69), as anticipated in the main text. 



More precisely, this equality has been derived for r) n -i = r) n = t]\ = 0. However, both the 



correlation function with Lagrangian insertions and the expression (88) are manifestly invariant 
under Z a -^x supersymmetries. Thus it follows that they agree for all values of the rj. 

E Relationships between different infrared regulators 

Several natural regulators are possible for the Wilson loop: 

• Dimensional regularization to D = 4 — 2e dimensions. 

• Dimensional regularization of the measure in the Lagrangian insertion procedure, keeping the 
integrand in four dimensions [27]. 



Multiplication of the Lagrangion insertion integrand by f] i ^j^wq^] (similar to adding 
masses in the Higgs-inspired scheme of 147] ) . 



and similarly several natural regulators are possible for scattering amplitudes: 
• Dimensional regularization to D = 4 + 2e dimensions. 

. MaMpeeatioe of the foot-di— integrand by Ul (sin* to above). 



Higgs regulator 47 



None of these regulators will produce the same answer. The last two differ in whether masses 
are kept in numerators. Here we would like to argue that the logarithm of the Wilson loop, which 
is the relevant physical quantity, can only differ from one approach to the other in a very mild and 
controlled way. 

More precisely, the logarithm of the Wilson loop takes the form (see for instance |5j and references 
therein) 

log(W n ) = £ ( 7 cus P log 2 + d log ™+I + C ) + Finite. (93) 

i 

In the scattering amplitude side this form follows from factorization theorems in the planar limit. 

It is completely expected, and verified in practice (to four loops in |48|), that different regulators 
only affect the scheme dependent constants d and C but the "finite part" is unchanged. This is part 
of the physical content of factorization. The modest point we would like to make in this Appendix 
is that factorization is made particularly explicit by the Lagrangian insertion procedure. 

First we observe that all the above regulators are compatible with the Lagrangian procedure (in 
fact, some are only defined within it). So, we can use it in our analysis. The main observation then 
is that in this procedure the logarithm can be taken already at the level of the integrand. That is, 
for instance, at two-loops, one can consider the integral 

log(W n ) = [ (M^( yi ,y 2 ) - ^Mr V (yi)M MHV ( y2 )) . (94) 
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The claim is that the parenthesis is much better behaved then each term in it object, with divergent 
regions with y\ and 1/2 not close canceling explicitly. This can be easily verified at two-loop using the 
well-known two-loop integrand for four particles |49|. More generally, we would like to argue that 
divergences in the logarithmic occur only from the regions where all particles collectively approach 
the same singular region. In other words, that factorization occurs already at the integrand level. 

The mechanism we suggest is very simple: Lagrangian insertions are color singlet operators. 
Imagine that, at t loops, a subgroup of Lagrangian insertions approaches a singular region. There is a 
scaling parameter s, for instance near a hard collinear singular region, {yf, y^ ,y^~) — > (yf, s 2 y~ , sy^~), 
which we can introduce to make all these insertions approach the hard collinear region at the same 
rate. Similar scaling variables can be introduced to deal with the soft wide-angle and soft-collinear 
regions. In the scaling limit we have a number of color-singlet operators inserted very close to an 
edge, plus other "passive" insertions at some other fixed locations. We know from general power- 
counting that the singularity in this limit will be at most logarithmic: the correlation function times 
the measure factor cannot not blow up in the limit. The only terms which survive thus behave like 
dimension-0 color-single corrections to the edge operator, as seen from the passive insertions. But 
such a correction can only be proportional to the edge itself. That is the field induced by the 
edge region is not affected, in the scaling limit, by the presence of the insertions. This implies the 
factorization of the integrand in this limit into a product of two correlation functions, one involving 
the singular subgroup, and one involving the passive subgroup. 

Since such products are killed by the logarithm, we conclude that the integrand for the logarithm 
of the W is regular whenever a proper subgroup of insertions approaches a singular region. As a 
corollary, if we take the scaling limit as all Lagrangian insertions approach a singular region, and 
strip off a collective d 4 y, the integral over the relative coordinates must be finite. This integrand 
multiplies the zeroth order expectation value of the Wilson loop. Since the mechanism is the color 
neutrality of the Lagrangian insertion, it applies in D dimensions as well and thus implies the 
equality of the first three regulators listed above, up to the finite "constants under the logarithm" . 
(The scaling limit in the soft directions, under which soft singularities factorize independently, 
is multiplicative. Thus the integral over relative rapidities, of the integrand for hard collinear 
divergences, is guaranteed to be finite.) 

A definite prediction of this argument is that there should exist a finite 4(1 — l)-dimensional 
integral which gives the cusp anomalous dimension to £ loops. Its integrand is the connected 
correlation function of I Lagrangian insertions in the universal background field very close to a 
segment of the Wilson loop, that is sourced by this segment and the two cusps at its endpoints. 
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